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I. INTRODUCTION 



Order and disorder parameters are useful tools in determining the phase structure of ex- 
tended physical systems. The order parameter of confinement in QCD is the temporal 
Wilson loop, which shows an area law in the confinement regime and a perimeter law in 
the deconfinement region [l| . The temporal Wilson loop is related to the potential between 
static quarks and an area law corresponds to a linearly rising quark potential. 

The emergence of the area and perimeter laws in the Wilson loop in the confined and 
deconfined phases can be easily understood in a center vortex picture of the Yang-Mills 
vacuum 2J. While confinement arises for a phase of percolated vortices, the perimeter law 
arises from small vortex loops, which can link to the Wilson loop only if they are in the 
vicinity of the loop. In recent years the center vortex picture of confinement has received 
strong support from lattice calculations [3|, |4j- 



Another (dis-)order parameter for Yang-Mills theory was introduced by 't Hooft [5j| and is 
defined in the following way. Consider Yang-Mills theory in canonical quantization (which 
assumes Weyl gauge A = 0), where the spatial components of the gauge field A(x) are 
the "coordinates". Let VKL4](C) denote the spatial Wilson loop operator. Then the oper- 
ator of 't Hooft's disorder parameter V(C) for a closed (spatial) loop C is defined by the 
commutation relation 

v(CiMc 2 ) = z L ^ c ^w{c 2 )v{c x ) , (i) 

where Z is a (non-trivial) center element of the gauge group and L(Ci,C 2 ) is the Gaussian 
linking number of the two loops C\ and C 2 . The expectation value of the (spatial) 't Hooft 
loop operator V(C) can serve as an order parameter for Yang-Mills theory. As argued by 
't Hooft, in the confined (deconfined) phase (V(C)} obeys a perimeter (area) law [6j. In this 
sense (V(C)) behaves dual to the temporal Wilson loop, which shows an area (perimeter) 
law in the confined (deconfined) phase. Given the fact, that a center vortex located at a 
loop C\ and described by a gauge potential A{C\) produces a Wilson loop C 2 

W[A{C X )]{C 2 ) = Z L{ - C ^ , (2) 
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the 't Hooft loop operator V(C) defined by equation ([!]) can be interpreted as a center vortex 
creation operator. An analogous monopole creation operator was studied on the lattice in 
ref. 3- 

In statistical physics, operators creating topological excitations like vortices or monopoles 
are referred to as disorder operators. Their expectation values, referred to as disorder 
parameters, are related to the free energy of the associated topological excitations. 

Center vortices can be easily generated on a lattice, where they represent co-closed (D — 2)- 
dimensional hyper-surfaces of plaquettes equal to a non-trivial center element. Accordingly, 
the 't Hooft loop operator can be easily realized on the lattice. In particular, center vortices 
wrapping around the whole lattice universe (torus), can be easily generated, and 't Hooft 
loops of maximal size on a finite lattice can be realized by imposing twisted boundary 
condition on the links [6]. The free energy of center vortices wrapping around the (space- 
time) torus, i.e. 't Hooft loops of maximal size, has been calculated in a high-temperature 
and measured on the lattice for SU(2) Yang-Mills theory in ref. 9| and for 



expansion 



a Z(3) random center vortex model [10(] in ref. Let us also mention that the energy 



per length of a straight infinite magnetic vortex line was calculated to one loop order as 



a function of the vortex flux 



It was found that the center vortex flux is energetically 



favored and its (free) energy vanishes in the absence of quarks. 

In the present paper we calculate the (spatial) 't Hooft loop, i.e. the vacuum expectation 
value of the operator V(C) defined by eq. (pQ), in continuum Yang-Mills theory. Recently 
progress has been made in variationally solving the Yang-Mills Schrodinger equation in 



Coulomb gauge [13| , 



14, LL 



161 ] . The Coulomb gauge has the advantage that Gauss' law 
can be explicitly resolved resulting in an explicit expression for the potential between static 
(color) charges. Using Gaussian ansatze for the vacuum wave functional, minimizing the 
vacuum energy density yields a set of Dyson-Schwinger equations (DSEs). These equations 



can be solved analytically in the infrared 17|. Implementing the so-called horizon condition 
(i.e. an infrared diverging ghost form factor) one finds a solution which produces a strictly 
linearly rising confinement potential for static color charges, provided the curvature of the 
space of gauge orbits, induced by the Faddeev- Popov determinant, is properly included. In 
this paper we use the Yang-Mills vacuum wave functional determined in refs. to 
calculate the spatial 't Hooft loop. 



The organization of the paper is as follows: In the next section we present the 't Hooft loop 
in Coulomb gauge, and the loop geometry will be worked out for a planar circular loo p. We 
will briefly review the Dyson-Schwinger equations in Coulomb gauge derived in ref. 14| in 
sect. IHII In sect. IIVI we study the asymptotic behavior for the 't Hooft loop in both the 
ultraviolet and infrared and extract its large perimeter behavior. Our numerical results will 
be presented in sect. |V] A short summary and some concluding remarks are given in sect. 

ED 



II. THE SPATIAL 'T HOOFT LOOP IN COULOMB GAUGE 
A. The 't Hooft loop operator 

The 't Hooft loop operator, V(C), is implicitly defined by eq. ([1]). An explicit realization of 
't Hooft's loop operator V(C) in continuum Yang-Mills theory was derived in ref. j3] and 
is given by 



V(C) = exp 



(3) 



where II" (x) = — i5/5A^(x) is the canonical momentum operator and 

^(S; x) = e J d 2 di5\x - x{a)) (4) 

is the gauge potential of a center vortex, whose flux is localized at the loop C = <9X. Here 
£ = i a T a (with T a being the generators of the gauge group) denotes a co-weight vector 
defined by exp(— £) = Z with Z being a (non-trivial) center element of the gauge group. 
Furthermore Xi(cr) denotes a parameterization of the 2-dimensional surface S bounded by 
the loop C. For SU (2) with generators T a = — |r a (r a being the Pauli matrices) the co- weight 
vector of the non-trivial center element Z = — 1 is given by £ a = 5 a3 2ir. 
Just as the Wilson loop W{C-i) creates an elementary electric flux along the loop C*2, the 
't Hooft loop V{Ci) creates an elementary magnetic flux along G\. In fact, the gauge 
potential (j3J produces a magnetic flux 

B(x) = d x A(E,x) =£ j dy5 3 (x-y) , (5) 

9S 
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which is localized on C = dH. 



Inserting (jlj) into ([2]), the spatial integral can be carried out yielding 



V(C) = exp 



(6) 



Since 11° (x) is the operator of the electric field, it is seen that the 't Hooft loop measures 
the electric flux through the surface £ enclosed by the loop C = <9£, in the same way as the 
(spatial) Wilson loop W(C) measures the magnetic flux through a surface enclosed by the 
loop C. 



The loop operator ([3]) is not manifestly gauge invariant. However, it produces a gauge 
invariant result, when it acts on a gauge invariant wave functional, see ref. 18] for more 
details. 



Consider the action of the loop operator |3| on a gauge invariant wave functional ^[A]. The 
't Hooft loop operator is analogous to the translation operator e m)3 in quantum mechanics 
which shifts the argument of a wave function exp(iap)if(x) = (p{x + a). Thus, when acting 
on wave functionals the loop operator ([3]) shifts the "coordinate" A(x) by the center vortex 
field A(x) 



V{C = dE)*[A] = V[A + A(Z)} 



(7) 



and in this sense is a true center vortex generator. 



Below, we firstly derive the expression for the 't Hooft loop in Coulomb gauge and secon dly 
calculate the 't Hooft loop for the specific vacuum wave functional obtained in refs. Q, 
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16l | . Furthermore we will explicitly work out the 't Hooft loop for a planar circular loop. 



B. The vacuum expectation value of the 't Hooft loop operator 

We are interested in the vacuum expectation value of the 't Hooft loop operator ([3]). The 
exact vacuum wave functional is not known but recently progress has been made in de- 
termining the Yang-Mills vacuum wave functional ^o(A) variationally in Coulomb gauge 
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OA = 



3, 



14 Il5l . Il6| . The wave functionals in Coulomb gauge satisfy Gauss' law 1 and 
hence should be regarded as the gauge invariant wave functionals restricted to transverse 
gauge fields. Using eq. ([7]) and implementing the Coulomb gauge by means of the Faddeev- 
Popov method, we obtain for the expectation value of the 't Hooft loop operator in Coulomb 
gauge wave functionals 

= J DA ± J(A ± )^*(A ± )^(A ± + A ± (Z)) , (8) 

where 

J{A L ) = Det{-Dd) (9) 

is the Faddeev-Popov determinant with D being the covariant derivative in the adjoint 
representation. 2 

Due to the transversal property of the field variable Af{x) only the transversal part 
A ± (dT,; x) of the center vortex gauge potential (TJl) enters, which is given by |20J 

Ai{dE, x) = -f j) da ik d%D(x - x(a)) , da ik = e iH dxi (10) 

ST, 

where D(x) is the Green's function of the 3-dimensional Laplacian defined by 

- A x D(x) = 6 3 (x) . (11) 
Note, that the transversal part A x) depends manifestly only on the boundary 



1 Gauss' law is explicitly resolved in deriving the Yang-Mills Hamiltonian in Coulomb gauge. 

2 The Fadeev-Popov determinant provides the iiaar measure of the gauge group. [l9( 

3 The longitudinal part of -4(£, x) is given by [20( 

A ll {E,x) = -£dn(E,x) , (12) 

where 

n(E,a?) = ( d 2 d k dtD{x-x(a)) (13) 



s 



is the solid angle subtended by the surface E at the point x. 



6 



f. Q 



In ref. [14| the following ansatz was used for the vacuum wave functional 

1 -Li / j3„, j3„/ A±a/ 



S^} = -Jd s xd 6 x'Af a (x)t i:j {x)u{x,x')Af a (x'), (14) 

where tij{x) = 5ij — dfdj/d% denotes the transversal projector and uj(x,x') is a variational 
kernel determined by minimizing the vacuum energy density. The ansatz (fT4l) for the wave 
functional is motivated by the form of the wave function of a point particle in a spherically 
symmetric (I = 0) s-state, which is of the form $?(r) = <ft(r)/r, where J = r 2 is the Jacobian 
corresponding to the change from Cartesian coordinates to spherically symmetric coordinates 
(For s-states the scalar product is given by (^1^2) — / drr 2 ^l(r)^2{i r ))- Like in the case 
of the point particle in a spherically symmetric state, the ansatz ffT4"j) simplifies the actual 
calculation, since it removes the Jacobian from the integration measure. Let us stress, 
however, that both, the ultraviolet and the infrared properties of the theory do not depend 
on the specific choice of the pre-exponential factor. In fact it was shown [l5[ , that replacing 
by J~ a with an arbitrary real a yields, after minimizing the energy density (to 2 loop 
level), a vacuum wave functional whose infrared limit is independent of the choice of a. (The 
choice a = | simplifies, however, the calculations.) 



For the wave functional ( 1T4|) the expectation value of the 't Hooft loop pj) becomes 

{V(dE)) = \Af\ 2 JvA ± J L 2[A ± }J-^[A ± + A ± }ex V (-(S [A ± ] + S [A ± + A ± })) .(15) 
15] it was shown, that to 2-loop level in the energy the Jacobian can be expressed as 
J (A) = const • exp f- j d 3 xd 3 x' A^^Uj^x^x, x^Af 11 ^')^ , (16) 



In ref. 



where Q 

5 2 In J{A~ 



1 1 sab. 



X^,y) = --— 2 — T 6 ab t kl (x)(t> 



4iV£-l 



¥) (17) 



5Ai«{x)5Al\y) 

is the "curvature" of the space of gauge orbits. This quantity gives the ghost loop contribu- 
tion to the gluon self-energy. The ghost propagator is defined by 

{^{-bdr^x'm = ^\^) (is) 

where d denotes the ghost form factor. 



Using the representation (|T6l) the functional integral in eq. (TT5T) becomes Gaussian and we 
obtain 



(V(dE)) = exp 



exp(-S) , (19) 



where we have introduced the abbreviation 



K(x,y)=u(x,y)-x(x,y)-~ / d 3 zdV(u;(x, z) - ^s, z 7 ) (a;(z', j/) - y))(20) 



If one uses the more general ansatz for the vacuum wave functional considered in ref. 151 ]. 

V(A) = J(A)- a e-^ AujA , (21) 

and uses the representation (fl~6l) for the Jacobian, one still finds the form (jl~9l) for the 't 
Hooft loop, however, with the kernel 

K = n-x-l(n-x)n-\n- x ), (22) 

where 

tt = uj + (1 - 2a)x (23) 

is the inverse (3-dimensional) gluon propagator. Note, eq. f l22|) arises from eq. fTSOj) with 
a; replaced by (|23|) . In ref. [15J] it was shown, that to one loop level (i.e. 2 loops in the 
energy density) the solutions of the Schwinger-Dyson equations for Q and x are independent 
of the choice of a. Thus, we obtain the same expectation value for the 't Hooft loop for all 
choices of a and for convenience we will continue to choose a = \ for which Q = uj. We 
should, however, stress that there are other quantities (like e.g. the 3-gluon vertex), which 
are sensitive to the choice of a fl7 |. 

Using eq. ( ITOi) one finds after straightforward evaluation for the exponent in the 't Hooft 

loop dm 

S = 2vr 2 j y k j dz x J d 3 xt kl (y)K(y, x)D{x, z) . (24) 
c c 

Due to the translational invariance of the vacuum the kernel K(x,y), as well as the Green 
function D(x,y), depends only on the distance \x — y\ and it is convenient to use Fourier 
representation 

D(x - x>) = J -0^ e ^-^D(q) , D(q) = ^ , t kl (q) = 5 kl - q kqi , q = j| . (25) 



Then the kernel ( 120]) becomes 

K{q) = (u(q) - X (q)) 
We obtain then from eq. 



x _ i ^(g) -x(g) 

2 u(q) 



oo 

S(C) = I dqK(q)H(C iq ) 



where the function 



H(C;q) = —^dy k ^ dz x / d^M^'^ 



(26) 



(27) 



(28) 



c c 

contains the geometry of the loop C considered, but is independent of the properties of the 
Yang-Mills vacuum, which are exclusively contained in the kernel K. 



C. Loop Geometry 

We can explicitly carry out the integral over the solid angle J dCl q in (j28~|h Expressing the 
unit vector q in terms of spherical coordinates G, 

q = q(Q, 4>) = (sin 9 cos 0, sin 6 sin 0, cos 6) (29) 

we have 



2tt 



r I 1 - ±sin 2 9 , k = 1,2 

/ d#«($(e,0 =27^ 2 (30) 

i I 1 - cos 2 9 , fc = 3 



Putting the 3-axis of g-space parallel to the vector y — z we find for the relevant integral 
n } fl( 1 + ^ 2 ) i jo(ot)-fo(a) , A; = 1,2 



^H(g)e mcosH = 27r»5 H / dze iaz { 2V = 2ir6 kl { (31) 



-i 



where a = q\y — z\ and we have introduced the spherical Bessel function jo(at) = 

To proceed further with the evaluation of the geometric function H(C, q) (128"]) we have to 
adopt an explicit realization of the loop C. For simplicity we will choose a planar circular 
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loop of radius R, which we put for convenience in the 1-2-plane. Then it is convenient to 
use cylindrical coordinates in x-space and to parameterize the loop by the azimuthal angle 

f 

x((p) = Re p (ip) , e p (ip) = ei cosy? + e 2 smip , dx = Re v {ip)dip , e v = —d v e p ((p) . (32) 
With eq. (13*11) the geometric function (1281) then becomes 

2tt 2tt 

H{C,q) = l -R 2 j dp J d^e^)-e^')f(Rq\e p (^-e p (^)\) , (33) 
o o 

where we have introduced the abbreviation 
i i 
f(a) = l -J dz(l + z 2 )e taz = l -J dz{\ + z 2 ) cos(az) = j Q (a) - j%(a) . (34) 
-l -l 



Using 



fo( x ) = -h( x ) , 3fi(z) = Jo(x) - 2j 2 (x) (35) 



the function f(x) fl34l) can be expressed as 

f( x ) = 2 -(2j,{x)- ]2 {x)) . (36) 

Since jo(0) = 1 and j«>o(0) = we have /(0) = | and the integrand in the angular integral 
is regular everywhere. Note also, that jo(x), j'o(x) and f(x) (IMl) are even functions of x. 

The expression ( l33i) can be further simplified by using 

e v (<p) ■ e<p(<P r ) = c os(<y3 - (p') , \e p ((f) - e p {(p')\ = a/2 (1 - cos(^ - <//)) (37) 
resulting in 

2tt 2tt 

H(C, q) = ^R 2 J dip J dip' cos(y? - (ff)f (Rqy/2(1 - cos(y? - tff) . (38) 



Hence the integrand in (|33l) depends only on the difference <p — (p'. Therefore we expect 
that one of the angular integrals can be trivially taken. This, in fact, turns out to be the 
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case, although the integrations run over finite intervals. In Appendix [A] we show that the 
geometric factor fl33|) can indeed be reduced to 



2 

H(C,q) = AnR 2 A(Rq) , A(x) = J da(l - 2 sin 2 a)f(2x sin a) , (39) 

o 

Let us stress, that the function H(C, q) reflects only the geometry of the considered loop C 
and does not contain any dynamical information as is clear from its definition given in eq. 
(128p . It contains only the kinemetical information about the transversality of the gauge field 
due the presence of the transversal projector. Note also, that the integrand in eq. (1391) is 
well-defined in the whole integration interval, in particular, the function f(x) is well-defined 
at the origin. After a sequence of manipulations (see Appendix [A} this integral (1391) can be 
done analytically with the result (x = qR) 

A(x) = — \2xJ (2x) + Tix (J 1 (2x)H (2x) - J (2x)H 1 (2x)) - 2J l (2x)} 
Ax 

[J 2 (2x)H 1 (2x) - J x {2x)H 2 (2x)\ , (40) 

where J u (x) and H u (x) are the ordinary Bessel functions and the Struve function, respec- 
tively. The asymptotic forms of the Bessel and Struve functions for small and large argu- 
ments are known. This yields the following asymptotic behaviors of the angular integral 
A(x) for small x 

A(x) = ^-x 2 - ^x 4 + • • • , x -> (41) 

and for large x 



A(x) — - — I- ^— cos ( 2x + — ) — = + O ( — = J 



x ->• oo. (42) 



Fig. [T] shows the full function A(x) together with its asymptotic forms. The asymptotic 
forms provide excellent approximations in the small and large momentum region. In the 
intermediate momentum range 1.5 < x < 2.5 we Taylor expand A(x) around x = 2 up to 
fourth order. Matching the Taylor expansion at x — 1.5 and x = 2.5, respectively, with 
the infrared and ultraviolet asymptotic forms yields a very accurate representation of the 
function A(x) in the whole momentum range, see fig. [TJ This representation is used in the 
numerical calculations. 
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FIG. 1: The angular function A(x) ()40f) : the full numerical result (solid line), its infrared (dashed 
line) and ultraviolet (short dashed line) asymptotic behaviors defined by eqs. (|4ip and (|42p and its 
Taylor expansion around x = 2 (dotted line). See text for more details. 

Inserting eq. (1591 into eq. (1271) we find for the exponent of the 't Hooft loop 

oo 

S(R) = AnR 2 J dqK{q)A{Rq) . (43) 
o 

We are interested in the large R behavior of S(R). For this purpose we study first the 
UV- and IR-behavior of the integrand, in particular of the kernel K(q). This requires the 
asymptotic behaviors of the solutions of the Dyson- Schwinger equations. 



III. THE DYSON-SCHWINGER EQUATIONS REVISITED 



A. Ultraviolet behavior and renormalization 



Consider the (unrenormalized) gap equation obtained in ref. 14 1 by variation of the energy 
density with respect to the kernel ui 

u;\k)=k 2 + X 2 (k)+Uk) + I°, (44) 
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where 

M-W'tJ^-W)*^ (45) 

^ =x/i^( i+ <^> 2 ) • g " - • - " x(q) y ' " — <46) 

are all UV-divergent integrals, so that this equation needs regularization and renormaliza- 
tion. In the above integrals d(k) and f(k) are, respectively, the ghost (see eq. ([IB]) ) and the 
Coulomb form factor, which are determined by two further DSEs, whose renormalization is 



described in detail in ref. 14] 



We renormalize the gap equation (1441) in the standard way by subtracting it at an arbitrary 
renormalization point fi. Defining 

*(*) = X(k) - = I x (k) - IM (48) 
and using xil 1 ) = we obtain 

u 2 (k) = k 2 + f(k) + oo 2 (v) - /i 2 + 2x(k) X (ji) + Uk) - I u (ji) . (49) 
The renormalized curvature x(k) ( )48l) is UV-finite. The crucial observation now is that in 



Iu>{k) P6l) the UV-divergent piece x(aO drops out. We can hence replace x{k) in Iu>(k) f|46j) 
by x(k). Analogously to ref. [ijj] we rewrite the Coulomb integral I u (k) (Hoi) as 

4W = ^ ) [x]W + 2x(A:)/«[x]W, (50) 

where the integrals 

N c f A I, , n ~>2\ d(k - q) 2 f(k - q) [uj(q) - X («)f - - *(*)]" 



«« = f/(^( 1 + ^ 2 ) 



(k-q) 2 w(q) 

(51) 

are for n = 1 and n = 2, respectively, linearly and quadratically ultraviolet divergent, while 
the differences 

A/j n) (*o = mm) - i^ixw (52) 

are ultraviolet finite. From eq. (1501) we obtain with x(a*) = 

/«(*) " W = UPik) + 2 x (k)l£\k) . (53) 
13 



With this result the renormalized gap equation (1491) can be rewritten as 

co 2 (k) - f(k) = k 2 + £ + Atf\k) + 2x(k) (e + AiJ)(fc) - A/W(0)) , (54) 

where we have introduced the abbreviations 

eo = ^ 2 (/i)-/i 2 (55) 
£ = x(A*)+^(0). (56) 

Since cu(/i) and /i are finite constants, £o is a finite renormalization constant. The only 
singular piece in the gap equation fl54|) is the quantity £ fl56l) . Ignoring the anomalous di- 
mensions both xif^) an d -fi^(O) are linearly UV-divergent. These linearly diverging terms of 
the 3-dimensional canonical Hamilton approach correspond to quadratically divergent terms 
in the 4-dimensional covariant Lagrange formulation. As is well known these singularities 
have to cancel by gauge invariance. We should hence ignore these singularities (which are 
artifacts due to our approximations) and keep only the finite parts of x(p) an d lu (0). In the 
following we therefore consider x(/ i )> and thus £ as finite renormalization constants. 

The novel result of the above consideration is that in the gap equation there are only two 
independent renormalization constants, £o and £. The renormalized gap equation obtained 



in ref. [14| seems to depend explicitly on the three renormalization constants: £, £o an d 
x(/i). The constant x(aO sneaked formally into the renormalized gap equation since it was 
not realized that this quantity drops out from the Coulomb integral I u (k) (l4"6l) . However, 
it was already found empirically in ref. [jjl that the numerical solutions to the DSE are 
practically independent of xilA- Indeed, by noticing that the integrals I^\x\{k) defined by 
eq. (15TT) satisfy the relations 

IL 1] [X + X^)}(k) = mm) (57) 
IL 2) [X + *(/*)](*) = Ii 2) [x](k) ~ MriWim) (58) 

one shows that the renormalized gap equation (1541) agrees with the one derived in ref. fl^ . 
Thus all results derived in ref. [l4] about the behavior of the solutions of the Dyson-Schwinger 
equations remain true. This refers in particular to the fact that both, the infrared and the 
ultraviolet behavior of the solutions of the gap equation are insensitive to the precise value 
of the renormalization constants (except for the renormalization constant of the ghost form 
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factor which is fixed, however, by the horizon condition, see below). While the (renormalized) 
gap equation obtained by minimizing the energy density does not depend on the choice of 
x(aO (but merely on the values of £ and £ ) other observables may explicitly depend on x (/•*)■ 
In fact, we will see that the 't Hooft loop does depend on x(a*)- 

The DSE resulting from the minimization of the energy density can be solved analytically 
not only in the ultraviolet where perturbation theory is applicable but also in the infrared 

3,0. 



B. Infrared behavior 



f. Q 



In ref. [141 it was shown, that the infrared behavior of the solutions of the DSE is rather 



uniquely determined once the so-called "horizon condition", [2l|], i.e. an infrared diver- 
gent ghost form factor, d~ 1 (k = 0) = 0, is implemented. This condition turns out to be 
absolutely necessary in D = 2 + 1-dimensions to obtain a self-consistent solution to the 
Dyson-Schwinger equation in Coulomb gauge [22]. In D = 3 + 1 this condition is in accord 
with the Gribov-Zwanziger confinement scenario, which is consistent with the Kugo-Ojima 



confinement criteria in Landau gauge 



23]. 



Infrared analysis of the ghost (or gluon) Dyson-Schwinger equation (in the rainbow-ladder 
approximation), implementing the horizon condition d~ l (k = 0) = and using the power 
law ansatze 

<"(*) = £ > « = p w 



yields the sum rule (in D = d + 1 dimensions) [Ti 

a = 2(3 + 2 - d (60) 

due to the non-renormalization of the ghost-gluon vertex. This sum rule also guarantees 
that x{k) (J35D (and thus x(&) dUD) has the same infrared power as u(k), i.e. 

u 

X(k) = -,k^0. (61) 

The horizon condition implies (3 > and for (3 > the gluon energy uj{k) (and thus also 
x(k)) is infrared divergent, i.e. a > 0. It can be shown that the gap equation has indeed a 
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consistent solution where both u(k) and are infrared divergent and in addition their 
difference (oo(k) — x(A;)) ^ s infrared finite. 

Indeed if (oo(k) —x(k)) is infrared finite it follows that l!?\k) and thus Al!?\k) are infrared 
finite. Then for infrared divergent u>(k) and x(/c) the gap equation (|54p reduces for k — ► 
to 



(u(k) + x(k)) (u(k) - x(k)) = ■ 

Since u;(fc) and x(fc) have the same infrared exponent it follows from 
limit of the gap equation is given by 

Jim (u(k) - x(k)) = £ 

or expressed in terms of the total curvature x(k) = X.(k) + xil 1 ) 

}im(u(k) - x(k)) = c , c = ^ - x(p) 

k— >0 



(62) 

that the infrared 
(63) 

(64) 



Instead of x(/x) we can use c as independent renormalization constant. Equation (|63|) implies 
(c.f. eqs. (1591) and (jBIl)) 



u = a 



(65) 



Equation (!63j) along with the ghost DSE in the infrared limit can be solved analytically for 
the infrared exponents yielding two solutions 

P w 0.796 , (3 = 1. (66) 

Only the latter one gives rise to a strictly linear rising static quark potential at r — > oo, 



V(r) 



2tt 2 



dk(d(k)f 



1 - 



sin kr 

kr 



(67) 



Here we have assumed that the so-called Coulomb form factor can be put to one [l4|. The 
first term in the bracket represents the self-energies of the static (quark and antiquark) color 
charges. In this paper we will use the solution (3 = 1 because the resulting linear rise of the 
static potential allows us to fix the scale by the (Coulomb) string tension. 



By definition the (Coulomb) string tension is given by that part of 

oo 

1 f „,,,/ (i sin x " 



dV 
dr 



2tt 2 



dkk(d(k)Y 



dQ(j 'tjc 



(68) 



x=kr 
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which (for large r) is independent of r. This part is exclusively determined by the infrared 
behavior of d(k). Inserting the infrared form fl59|) (with (3 = 1) into eq. fl68l) we obtain for 
the Coulomb string tension 



This integral can be done, yielding 



2tt 2 



dx d ( sinx 

x dx 



x 



(69) 



57T 



(70) 



Thus the infrared coefficient of the ghost form factor can be explicitly expressed by the 
string tension. The infrared analysis of the ghost Dyson-Schwinger eqnation provides the 
relation between the infrared coefficients in fl5T?|) [17] 



b 2 N c r(^) 2 r(i-f + p) 

4(4tt)^ T(d-P)T(l + §) 2 



(71) 



Hence the infrared coefficient a of u(k) and x(k) (see eqs. ( |59l) . ( ISTl) . ( 1651) ) can be also 
expressed by the string tension. For d = 3 and (3 = 1 the last relation simplifies to 



8tt 



2 ' 



(72) 



Lattice calculations 



241 ]. 251 ] show that o"c is by a factor of 1.5 - 3 larger than the string 



tension a extracted from the Wilson loop, in accord with the fact that ac gives an upper 
bound to cr 26]. Assuming a value of the Coulomb string tension ac = 1.5<t(ct = (440MeV^) 2 ) 



we find 



a = 0.185 GeV 2 



b = 2.702 GeV 



(73) 



These values are in remarkable agreement with the figures extracted from the numerical 
solutions of the Dyson-Schwinger equations. 

IV. THE ASYMPTOTIC BEHAVIOR OF THE 'T HOOFT LOOP 



Given that uj(k) and x(/c) are both infrared singular and satisfy the relation (|64l) it follows 
that the infrared limit of the kernel K(k) (T26]) is given by 



lim K(k) = lim(u;(A;) — x(^)) — c 



(74) 
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As we will explicitly see further below, it is the infrared behavior of K(k) which determines 
the large i?-behavior of the 't Hooft loop. The fact that K(q) is infrared finite (see eq. (1741) ) 
restricts the possible independence of S (H5|) drastically. 



A. Renormalization of the 't Hooft loop 

In ref. [ijj it was shown that in the ultraviolet (k — > oo) 

v(k) - k , - -4= > ( 7s ) 

where /i is some arbitrary energy scale (renormalization point) and c' is some constant whose 
precise value depends on the choice of \x. The first relation expresses asymptotic freedom, i.e. 
for large momenta gluons become photons, while the second relation expresses the presence 
of an anomalous dimension. With ( 1751) for large q — > oo the kernel K(q) defined by eq. 
( 125]) behaves asymptotically as K(q) — > \q. Since A(x = qR) ~ - for x — > oo, see eq. 
(|42|) . it follows that the exponent of the 't Hooft loop S (jl5|) is linearly UV- divergent. UV- 
singularities are expected given the singular nature of the 't Hooft loop. There are also 
UV-subleading terms in K(q) which give rise to UV-singularities in S (H3|) . To extract these 
terms let us parameterize ui(q) by 

u{k) = ^ + a + k, (76) 

which is in accord with the asymptotic IR- and UV-behavior and which gives a good approx- 
imation to the numerical solution of the gap equation in the whole momentum range (see 
fig. [2]). For the solution presented in fig. [2] we find numerically a = 0.61 and clq = —0.078. 

The leading UV-terms of K(q) giving rise to UV singularities in S (H3l) are then given by 



= \ ( i-\^ + \^ + c^)-a ^ I + ... . (77) 

with unknown coefficicents c, c_i arising from the UV-expansion of x(l)- These terms 
produce UV-singularities of S ( 1431) of the form (A - UV-cut-off) 

~ A, li(A) , In A, ln(lnA) , (78) 
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r^IG. 2: The gluon energy resulting from the numerical solution of the Dyson- Schwinger equation 
1^ ] and its parameterized form (|76|) . 

where \i(x) is the logarithmic integral, which behaves for large x as li(x) = tt(x)+0 (y/x lnx), 
where ir(x) is the number of primes less than or equal to x. In practice it is not possible 
to determine the various coefficients of the UV-asymptotic expansion of x{q) with sufficient 
accuracy. We therefore propose to eliminate the UV-divergencies of S ( 1431) by subtracting 
from the kernel K(q) ( 1261) the following UV-leading one 

A'ofe) = ife + O «)(l-(^|) 2 ) , (79) 

which results from K(q) ( 1261) by replacing u(q) by its UV-leading and next to leading order 
terms (see eq. fl76|) ). while keeping the ratio x(q) /w(q) fixed. The kernel K (q) fl79l) contains 
all the UV-leading terms which give rise to UV-singularities in S. In addition (since we 
have kept fully ^1 — (J) 2 )) it also contains terms which produce finite contributions to S. 
However, these terms are negligible compared to the dominating finite terms of S produced 
by the full K. The reason is the following: For large R the 't Hooft loop is dominated 
by the infrared part of K(q) as one explicitly notices by rescaling the integration variable 
q — > qR = x in (H3|) and as we will explicitly show further below. Furthermore, since 
u(q) ~ 1/q for q — > 0, in the infrared K (q) is by a power of q suppressed relative to 
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K(q) and thus gives sub-leading contributions to the large R behavior of the 't Hooft loop. 
By the same token, elimination of the UV-divergent terms from S by replacing K(q) by 
K(q) — K (q) does not change the large i?-behavior of S. 

We thus isolate the ultraviolet divergent parts of S (H31 by writing 

S = (S - S ) + S , (80) 

where So results from S (H3l) by replacing K(q) (|26|) by K (q) (1791) . By construction S — So is 
UV-fmite. The leading UV-divergent term S^ of So is independent of R and exp(— Sq ) = Z 
can be absorbed into the renormalization of the wave functional Z^ty — > The remaining 
terms in Sq have weak independences (at most logarithmic) and as explained above their 
elimination does not spoil the use of the renormalized 't Hooft loop 

oo 

S = S-So = AkR j dx(K (|) - Ko (I)) A{x) (81) 
o 

as a disorder parameter, i.e. they do not convert a perimeter law into an area law or vice 
versa. Eq. ( 1HT1) explicitly shows that the large i?-behavior of the 't Hooft loop is determined 
by the infrared behavior of K(q). 



B. The large R behavior of the 't Hooft loop 



We are interested in the 't Hooft loop, i.e. in the quantity S(R) defined by eq. (TSTT) for large 
R — > oo. If we could take in eq. ( IHTl) the limit R — > oo before taking the momentum integral, 
using 



lim K (^) = K(0) = c 
R^oo \RJ 

we would obtain a perimeter law 

oo 

S = Rn, k = 4tic / dxA(x) (83) 



o 



for c 0. However, the remaining integral (1831) is UV-divergent while the original integral 
f lHT]) is convergent. Thus we must not take the limit R — > oo before having carried out the 
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momentum integral: In fact, we will show below that a finite c ^ gives rise to a large 
i?-behavior S(R) ~ R\nR, while c = yields a perimeter law S(R) ~ -R. 



The fact that we must not take the limit i? — > oo in (IHTj) before taking the integral does 
not, however, imply that S 1 is not dominated by the infrared part of K(q). We will explicitly 
show further below that it is the infrared part of the momentum integral (l8Tj) (and thus of 
K(q)) which determines whether an area law or perimeter law arises. 

Since u{k) — x{fy is infrared finite (by the gap equation, see eq. ( 1631) ) it is of the form 

uj(k) - X (k) = CkT 1 , 7 > , k -»• , (84) 

and by eq. (|74p . K(k) has the same IR-behavior 4 

K{k) = CW , 7 > , k -»• . (85) 

Since the large i?-behavior of S(R) is determined by the infrared behavior of K(q) and since 
the UV-behavior of K(q) and A(Rq) determines the convergence property of the momentum 
integral, for the following qualitative considerations it is sufficient to approximate A(x) by 
the interpolating formula 

a, s 27T x 2 8 

15 r/x 6 + 1 15 

which has the correct IR and UV behavior given in eqs. (T4T!) and (|42|) . respectively. Along 
the same lines we use the interpolating formula 

AK (,) := K(„) - K ( 4 ) = (87) 

which has the correct IR and UV behavior. With eqs. flH7|) and flHBT) we find from ([8*1]) 



oo 



-> S = / rf s -v+i I 77 pwu I 1 = : J ^ + T W ■ (88) 



8?r 2 Ci2 7 x^ +1 + (g ^) 7+1 ^/^ 3 + 1 



o 



To estimate this integral we split the integration range into the intervals [0, q R] and [qoR, oo) 
and denote the corresponding integrals by Ijr and Iyy. Consider first Ijr where < x < 



In the notation of eq. (|74|) 7 > corresponds to c = and 7 = 0toc = C7^0. 
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q R. An upper and lower bound to this integral is obtained by replacing x -y+i + ( qoR ^+i by 



(go/?) 



7 + 1 



and 



(q R)->+ 1 +(qoR)~>+ 



i = I ( gofi )7+i ; respectively. Thus we obtain 



(89) 



where 



1 IR 



q R 



qoR 



dx- 



X 



7+2 



?7a; 3 + 1 



An upper and lower bound to the integral Ijjv is obtained by replacing 



z 7+1 +(<20^) 7+ 



(90) 



t by ^ 



and 



1 



Xf + 1 +XT + 1 2x~t + 



i+r, respectively. This yields 



(91) 



where 



l uv 



qoR / <ix 



x 



goR 



t]X 3 + 1 



3^7" 



V3 ( arctan - _ 

2 y/3 



1 y 2 + 2y + 1 
H — In — 5 

2 - y + 1 



(92) 



Here we have introduced the abbreviation y = i] 1 ^ 3 q R. Note that is independent of 7 
and thus independent of the IR-behavior of AK(q). For large R the last expression becomes 



1 uv 



1 

n 



1 - 







R 5 



(93) 



4?7 (q R) 3 _ 

Thus the UV-part contributes for large R always a perimeter term to 5 flHHl) . It is the 
IR-part, or more precisely the infrared exponent 7 of K(q), which decides whether the 
't Hooft loop has an area law or a perimeter law. 



Consider first 7 = 0. In this case the IR-integral (190]) becomes 

gR 



1 IR 



dx- 



xr 



1 



T = -ln [1 + 77(q ^) 3 ] 



rjx 3 + 1 3rj 



(94) 



and gives rise to an .R-dependence 



S ~ RlnR, R->oo. 



(95) 
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Consider now the case c = 0, i.e. 7 > 0. For 7 > the IR-integral (!90|) can be expressed as 



Iir(R) 



—IT 



37](rj 1/3 q Rysm — 



7T7 



-1 



377 



r]{q Rf 



i,- 5 7 



where 



$(z,s,a) = 



k=0 



a + k)' 



(96) 



(97) 



is the Lerch transcendent, which is a generalization of the Hurwitz zeta function and the 
poly logarithm function and is defined for \z\ < 1 and a 7^ 0, —1, ... . For large R we can 
expand &(z, 1, —§7) in powers of z = — l/(r]q R) 3 , keeping only the first two terms in ([9j 
In this order we obtain 



Iir(R) 



-71 



3r] [r] 1 q R) sin — 



1 1 1 

777 3r]7](q R) 3 (l 



h) 



(98) 



For 7^0 the first two terms have simple poles which cancel. For 7 > the second term is 
the leading one for large R, and gives rise to a perimeter contribution to S fl88l) . Since the 
UV integral (|93|) contributes a perimeter term independent of 7 we obtain for 7 > indeed 
a perimeter law for the 't Hooft loop. 

We have thus shown that the Yang-Mills vacuum wave functional obtained in the variational 



161 ] yields indeed 



solution of the Yang-Mills Schrodinger equation in Coulomb gauge 
a perimeter law for the 't Hooft loop, provided the renormalization constant (1711) is chosen 
c = 0. The 't Hooft loop is thus more sensitive to the details of the wave functional than 
the confinement properties, i.e. the infrared behaviors of the ghost and gluon propagators, 
which turned out to be independent of the choice of the renormalization constants once 
the horizon condition is implemented [if]]. The choice c = is also preferred by the 
variational principle as we will discuss in sect. IVI1 



V. NUMERICAL RESULTS 



The above observed qualitative behavior of the 't Hooft loop is confirmed by the numerical 
calculations. The DSEs resulting from the variational solution of the Yang-Mills Schrodinger 
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FIG. 3: The difference ui{k) — x(fc) for various choices of the renormalization constants. For all 
solutions presented here the renormalization constant c = (see eq. (I64p ) is chosen, resulting in 
an infrared- vanishing u(k) — x(k)- 



equation in Coulomb gauge derived in ref. 14J and reviewed in sect. IHII are solved numerically 
as described in ref. 16]. 

Since Yang-Mills theory is a scale-free theory, the DSEs can be entirely expressed in terms 
of dimensionless quantities denned by rescaling all dimensionful quantities with appropriate 
powers of a (so far arbitrary) scale \i: 

u)(k) = uj(kfi)/fi, x(k) = x(kfi)//J,, dik) = d(kfi), f(k) = f{kp). (99) 

From the numerical solution of the DSEs we calculate the dimensionless Coulomb string 
tension a c (169|) . which is related to the physical (Coulomb) string tension cr c by 



Or 



(100) 



The last relation fixes the scale \i once a specific value is assigned to a c . Note that a c 
depends, in principle, on the renormalization constants £, ^o- 

Fig. [3] shows the difference u(k) — x(k) for various choices of the remaining renormalization 
constants £, in table [J we quote the resulting dimensionless Coulomb string tension a c . 
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0.0 


0.0 


0.0 


1.44 


1.44 


1.44 




-0.15625 





0.15625 


-0.15625 





0.15625 




0.047 


0.173 


0.340 


0.261 


0.385 


0.285 



TABLE I: The dimensionless renormalization constants £, £o used in the numerical calculations 
presented in fig. [3J and the resulting dimensionless Coulomb string tension a c . 

For all these solutions the kernel K(q) (|26|) of the 't Hooft loop has practically the same 
infrared behavior (JT4"|) . fl84"|) with 7 « 1. With these solutions the quantity S'(-R) flHTj) of the 
't Hooft loop is calculated. Thereby we use for the angular integral A{x) (HOI) the analytic 
representation presented at the end of sect. [TT1 In the numerical evaluation of the momentum 
integral in S f lHTj) we split the UV part of A(x) (see eq. (1421) ) into an oscillating and non- 
oscillating part. The non-oscillating part as well as the infrared part is integrated by the 
standard Gauss-Legendre method. The integral over the oscillating part of A(x) is efficiently 
carried out using the method of integrating half-waves. Fig. 4 shows the results obtained 
for S(R)/R flSTl) for c = and various values of the remaining renormalization constants 
£, £o- One observes that indeed for large R the quantity S(R)/R approaches a constant 
independent of the choice of the remaining renormalization constants while requiring c = 0. 
Although all these solutions with c = yield a perimeter law for the 't Hooft loop, i.e. 
S(R) ~ R, the perimeter tension k defined by 

S(R -> 00) -> kR (101) 

depends quite sensitively on the remaining renormalization constants £, £o- m f &c t these 
renormalization constants affect the solutions of the DSE only in the intermediate momen- 
tum range while both the infrared and ultraviolet behaviors are untouched. This shows that 
the perimeter tension is sensitive to the details of the wave functional in the intermediate 
momentum range. Therefore it would be very helpful to have lattice measurements of this 
quantity. Unfortunately this quantity is very difficult to measure on the lattice. 

The range of the undetermined renormalization constant implies quite large changes of the 
gluon energy in the intermediate momentum range. If we interprete the lowest-lying glueball 
as a "two quasi-gluon s-state" and adjust the minimum in 2u(k) to the glueball energy of 
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1.65 GeV and furthermore choose the Coulomb string tension 1.5 o [o = (440MeV) 2 ) we 
find the perimeter tension around 4.4 GeV. 
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FIG. 4: The function S(R)/R, (|8ip . for different renormalization constants, resulting from the 
solutions presented in fig. [3l 



VI. CONCLUDING REMARKS 



In summary, the wave functional we have obtained in the variational solution of the Yang- 
Mills Schrodinger equation yields either a perimeter law 5 ~ R or S ~ i?ln R depending 
on whether we choose the renormalization constant ( 1741) c = or c ^ 0, respectively. The 
renormalization constants are undetermined parameters due to the approximations involved 

n 

in the derivation of the DSEs [14J]. However, since we are using the variational principle 
the renormalization constant could, in principle, be considered as variational parameters 
and fixed by minimizing the energy density with respect to these parameters. Such type of 
calculations are rather expensive, since they require the solution of the DSE in the multi- 
dimensional space of the renormalization constants and as the energy density depends only 
implicitly via the various form factors on these constants. An exception is the renormal- 
ization constant c defined in eq. fl64"|) . which explicitly enters the kinetic energy and which 
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determines the large i?-behavior of the 't Hooft loop. Since c is an infrared property of the 
solutions of the DSE for its determination it should be sufficient to use the infrared form 



of these solutions which are analytically known [17j. Then c drops out from the Coulomb 
energy and does not enter the potential energy, while the kinetic energy density is propor- 
tional to c 2 and is minimized by c = 0. It is this value which gives rise to a perimeter law 
in the 't Hooft loop. Resorting to the representation ([16]) (which is correct to two loop in 
the energy) the value c = implies that our wave functional (THl) has the form 

^(A)=Afe-y AhA (102) 

with 

h{k) = u(k) - x{k) ~ k 1 , k -»• , j > . (103) 

It describes a vacuum where the low momentum modes of the gauge field can fluctuate (for 
k — > 0) in an unconstrained fashion. Such a stochastic vacuum is required for an area law 
in the spatial Wilson loop. We therefore find that a perimeter law in the 't Hooft loop is a 
necessary condition for an area law in the spatial Wilson loop, which at zero temperature 
is equivalent to the temporal one. This is in accord with the duality between the spatial 
't Hooft loop and the temporal Wilson loop. 
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APPENDIX A: THE LOOP INTEGRAL 



In the following we work out the geometric factor H(C, q) given by ( 1381) of the 't Hooft loop. 
We first reduce the double integral to a single loop integral, which eventually can be taken 
explicitly. 

Since the integrand in eq. ( 1381) depends only on the difference ip — <p' it is convenient to 
introduce "center of mass" and "relative" coordinates 

ip + tp' 

(f) = — - — ,*y = <p-<p . (Al) 
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The integral (|38|) then becomes 
1 



H(C, q) = -R 2 



cos7/( J Rg v / 2( 1 - cos 7)) . (A2) 



n 20 2tt 2(2tt-0) 

-20 7T -2(2vr-0) 

Changing integration variable in the second term — * 7r + and using cos(— 7) = cos 7 we 
find 

20 2(tt-0) 



H(C,q) = R 2 



cos-yf(Rq^/2(l — 00S7)) 



(A3) 



c?7 + y c?7 
.0 

A further change of integration variable 7 — > 7 — 2ir in the second term using the periodicity 
of cos 7, yields 



H(C,q) = R 2 



- 20 


■ 20 -20 _ 




h 


J c?7 + / c?7 




.0 


.0 -2tt 





cos7/(i2gy2(i — C0S7)) 



7T 27T 



= R 2 J J d>ycosif(qRy/2(l - cos 7)) 


2tt 

= ixR 2 J d-f cos j f(qR^ 2(1 - cos 7)) . (A4) 


As anticipated, one angular integral can be trivially taken. Changing the integration variable 



a = ^ we obtain 



F(C, g) = AnR 2 A(Rq) 



(A5) 



where 



v4(x) = / da(l — 2 sin 2 a)f(2x sin a) 




(A6) 



which is the result quoted in eq. ([39 



To work out the angular integral (1A6jl . it is more convenient to express /(a) ( |34l) in terms 
of its original integral representation, rather than in terms of the spherical Bessel function, 

1 

(A7) 



/(*) = ~ y dz (1 + cos(te) . 



-1 
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Then, the angular integral in (1A6[) can be explicitly carried out. This turns out to be also 
advantageous for the numerical calculations. 

To carry out the angular integration in eq. (1A6j) we use in /(2xsina) (IMI) the expansion 

OO 

cos(y sin a) = Jo(y) + 2 J2n(y) cos 2na , (A8) 

n=l 

where J v (x) are the (ordinary) Bessel functions. 
Then we need the following angular integrals 



da cos 2na = — sin 2na 
2n 



for n 7^ 



o 



2 2 



c/(i2 sin' n cos2no! = J <ia;(2sin a — 1) cos2raa 
o o 



2 



7T 



<ia cos 2a cos 2na = — — 8 Ui i 



(A9) 



(A10) 



To obtain the last relation we have inserted a zero in form of the integral (1A9I) . With eq. 
flA9l we find from flA8|) 



7T 



da cos(?/ sin a) = —Jo(y) 



(All) 



Analogously, we find from eq. (1A8I) with (lAlOj) 

2 

y da2 sin 2 a cos(y sin a) = —(Jo(y) — 'h{y)) ■ (A12) 
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With the last two relations we obtain for the angular integral (1A6 



2 



A(x) — - f dz(l + z 2 ) J da(l — 2 sin 2 a) cos(2x;z sin a) 



-l o 
i 



= i J dz{\ + z 2 ) [^J (2xz) - | (J (2xz) - J 2 (2x2;)) 
-l 

i 

= \J dz(l + z 2 )J 2 (2xz) . (A13) 
-l 

The remaining integral over the Bessel function can be taken using the relations 

z z 

J dtJ 2 (t) = J dtJ (t)-2J 1 (z) 



z 

j dtJ (t) = zJo(z) + ^z(H (z)J 1 (z) - H 1 (z)J (z)) 



1 

/3% 
dzz 2 J 2 (zy) = — [J 2 {y)H x {y) - H*{y)J x {y)\ , (A14) 

o 

where H v (z) are the Struve functions. 

With these relations we find for the angular integral (1A6|) the result quoted in eq. (H0|) . 

We will also need the asymptotic form of A(x) for small and large arguments. This is most 
easily done by using the asymptotic forms of the Bessel function J 2 {x) in (IA13I) (although 
the asymptotic form of the Struve function is also known). One finds then the results quoted 
in eqs. (Hill and (jl2l). 
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